Fast two-bit operations in inductively coupled flux qubits 
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A central problem for implementing efficient quantum computing is how to realize fast operations 
(both one- and two-bit ones). However, this is difficult to achieve for a collection of qubits, especially 
for those separated far away, because the interbit coupling is usually much weaker than the intrabit 
coupling. Here we present an experimentally feasible method to effectively couple two flux qubits 
via a common inductance and treat both single and coupled flux qubits with more realistic models 
which include the loop inductance. The main advantage of our proposal is that a strong interbit 
coupling can be achieved using a small inductance, so that two-bit operations as fast as one-bit ones 
can be easily realized. We also show the flux dependence of the transitions between states for the 
coupled flux qubits. 

PACS numbers: 74.50.+r, 85.25.Cp, 03.67.Lx 



I. INTRODUCTION 

Josephson-junction circuits can exhibit quantum be- 
haviors. Among qubits based on Josephson-junction cir- 
cuits, the charge qubit realized in a Cooper-pair box can 
demonstrate quantum oscillations^ An improved version 
of this circuit has showed quantum oscillations with a 
high quality factori* In addition to charge qubits, flux 
qubits achieved in a superconducting loop with ona 3 - or 
three Josephson junctions* have been studied and some 
of these have shown quantum dynamics A The phase 
qubit consists of a large-area current-biased Josephson 
junction^ 

Capacitive couplings of two superconducting qubits 
(both charge— and phase-types*) were attained recently 
in experiments, and quantum entanglement was observed 
in these systems. Also, controllable interbit couplings of 
charge qubits were proposed using a variable electrostatic 
transformer j* a current-biased Josephson junction 10 and 
a tunable dc-SQUTDii* These interbit couplings can 
link nearest neighboring qubits. Actually, there are 
quantum-computing protocols (e.g., adiabatic quantum 
computing 1 -) that only demand nearest-neighbor cou- 
plings. However, for more general quantum-computing 
protocols, it is desirable to achieve strong enough cou- 
plings among non-neighboring qubits as well. When 
charge qubits are coupled by LC-oscillator modes* 3 - or 
by an inductance^* long-ranged interbit couplings can 
be realized, but a very large value of the inductance is 
needed. An alternative way of coupling charge qubits was 
proposed using a Josephson junction^i*** Moreover, 
the charge qubit can be very sensitive to the background 
charge fluctuations, which generate noise that severely 
limits the performance of charge-qubit devices and, un- 
fortunately, is difficult to reduce. 

In this paper, we present an experimentally feasible 



method to effectively couple two flux qubits. In con- 
trast with the charge qubit, the flux qubit is insensitive 
to the charge noise. In this qubit, the major noise is 
due to the fluctuations of the magnetic fluxes. Estima- 
tions show that the flux qubit can have a relatively high 
quality factorii* Here we include the effect of the loop 
inductance in a three-junction flux qubit and couple two 
flux qubits via a common inductance. Because the criti- 
cal current of each Josephson junction in the flux qubit 
is larger than that in the charge qubit, we can produce 
a strong interbit coupling using an inductance as small 
as 20 pH (corresponding to a loop diameter of approxi- 
mately 16 (im and comparable to the loop inductance of 
the single flux qubit currently achieved in experiments), 
and thereby two-bit operations as fast as one-bit ones can 
be easily achieved, improving the efficiency of quantum 
computing. Moreover, we show a novel flux dependence 
of the state transitions in two coupled flux qubits. We 
find that, except for some specific values of the external 
flux, the forbidden transitions in the two coupled flux 
qubits become allowed when the parameters of the two 
qubits change from being initially equal to each other and 
then making these different. 

Coupling two fl ux qubits by a mutual inductance was 
proposed in Refs. Il9f20|2ll and was recently realized in 



experiments. 



Here we treat both single and coupled 



flux qubits using more realistic models which include the 
loop inductance. We numerically solve the Schrodinger 
equation to obtain the energy levels and the eigenstates 
of the flux-qubit systems. This numerical method allows 
us to extend our study to the larger inductance regime. 

The paper is organized as follows. In Sec. II, we study a 
single flux qubit containing loop inductance. It is shown 
that the system can still be used to achieve a qubit even 
for a larger loop inductance of L ~ 1 nH. Section III 
focuses on two flux qubits coupled by a common indue- 
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FIG. 1: (Color online) (a) A flux qubit, where an exter- 
nal magnetic flux $ e pierces the superconducting loop that 
contains three Josephson junctions and an inductance L. 
The Josephson energies and capacitances of the junctions 
are Eji = Ej 2 = Ej, Ci = C 2 = C, E J3 = aEj, and 
Ci — aC. Here we choose a — 0.8 and E.j — 35E C , where 
E c — e 2 /2C. (b) Two flux qubits coupled by a common in- 
ductance L c , where the external flux $ e is applied within 
the left loop A\L c B\Ai. The parameters of each flux qubit 



are = E% = E, 
and C { i ] 
and Ef 



(i) 



= C* w , £ 



with i = 1,2. Here we choose ai = 0.8 
where ^ = e 2 /2C (l) . To implement a 
readout of the flux-qubit states, a switchable superconducting 
flux transformer is employed to couple the dc-SQUID magne- 
tometer with the inductance L in (a) or L c in (b) during the 
quantum measurement. However, this coupling is switched 
off in the absence of a readout. 



tance. In Sec. IV, we study the state transitions induced 
by the microwave field. Section V deals with the circu- 
lating supercurrents and quantum measurement. Finally, 
the discussion and conclusion are given in Sec. VI. 



II. SINGLE FLUX QUBIT 
A. The model 

We first consider a single flux qubit in the absence of 
a quantum measurement, where the dc-SQUID magne- 
tometer for measuring quantum states of the flux qubit 
is decoupled from the qubit. As shown in Fig. 1(a), 
the flux qubit consists of a superconducting loop with 
three Josephson junctions and the total inductance of 
the whole loop is L. Fluxoid quantization around the 



loop imposes a constraint on the phase drops across the 
three junctions: 
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Here, $o = h/2e is the flux quantum, 

/ = *e/*0 

represents the reduced magnetic flux, and 
I = losing, 



(1) 



(2) 



(3) 



(4) 



with Iq = 2ttEj /$0j is the circulating supercurrent. 

When the loop inductance is included, the Hamiltonian 
of the single flux qubit is 



p2 

H= p 



P 2 



2M P 2M„ 



U{<j> p ,(j> m ), 



(5) 



with the potential energy given by 

U(<f>p, 4> m ) = Ej[2 + a- 2 cos (j> p cos <f> m 

-acos(27r/' + 20 m )] + ^LI 2 . 



(6) 



Here 



Pk 



■ih- 



p,m, 



M p = 2C($o/27r) 2 , 
M m =M p {l + 2a), 



(7) 



and 
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4>i ~ 4>2)- 



Also, the supercurrent / can be rewritten as 
/ = I sm.(<p p +<f> m ). 



(8) 



(9) 



The Hamiltonian (5) is reduced to Eq. (12) in Ref. Il9l 
when L — > 0. 

Figure 2 presents the contour plots of the periodic po- 
tential U(</)p,(f> m ) for / = 0.5 and a = 0.8. The nu- 
merical results show that the minima of the potential 
preserve the two-dimensional centered cubic lattice even 
for a large loop inductance. For inductance ratio 



0l = L/Lj 



(10) 



from zero to one (where Lj = $ /2irI is the Josephson- 
junction inductance), a well defined double- well poten- 
tial structure exists at each lattice point even though at 
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FIG. 2: (Color online) Contour plots of the potential energy 
U(<f> p ,4>m), in units of Ej, for a = 0.8 and / = 0.5. Here 
f3 L = 27r7 L/$ = (a) 0, (b) 1, (c) 4, and (d) 10. Notice 
that the well defined double-well potential structure vanishes 
in (d), and thus the flux qubit breaks down. 
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FIG. 3: (Color online) Energy levels of a single flux qubit 
versus reduced flux / for different values of (3l, where only 
the levels of the states i = to 5, are shown. Here the 
energy E is in units of Ej. Notice the robustness of the two 
lowest levels for wide changes in the loop inductance L. 



higher energies the well shapes are modified by the loop 
inductance L. This double- well structure is required for 
achieving a two-level system. As shown in Fig. 3, the 
lowest two levels of the single-qubit system are not signif- 
icantly affected by the variation of (3l (when < Pl ^ 1) 
because the corresponding two eigenstates are mainly 
contributed by the weakly /^-dependent ground state in 
each well. However, since varying L significantly modi- 
fies the well shapes at higher energies, the excited states 
within or above the wells (which, as seen in Fig. 3, dom- 
inantly contribute to the eigenstates corresponding to 
the third and higher levels) become pronouncedly /Si- 
dependent. Indeed, Figure 3 shows that the top three 
levels are sensitive to the variation of (3l (even when 
^ Pl ^ 1)- Moreover, with the loop inductance in- 
creasing to (3l ~ 4 [see Fig. 2(c)], a more distorted 
double- well structure appears at each lattice point, and a 
local energy minimum develops along the diagonal direc- 
tion between every two adjoining double-well structures. 
These newly-developed local minima will affect the two- 
level system achieved for the qubit. When the loop in- 
ductance increases even more to (3l ~ 10 [see Fig. 2(d)], 
the periodic potential is even more distorted. In this case, 
the well defined double- well potential structure vanishes, 
and thus the flux qubit breaks down. 



B. Energy spectrum 

The energy spectrum and the eigenstates are deter- 
mined by 

H^{(/> p ,(/> m ) = E^((f )p ,(f )m ). (11) 

Figure 3 shows the dependence of the energy levels on the 
magnetic flux for fti, < 1. Here we choose Ej = 35E C , 
where the charging energy E c is defined as E c = e 2 /2C. 
These parameters are close to those used in a recently 
fabricated flux-qubit device £ 

Around / = 0.5, in sharp contrast with the higher 
energy levels, the energy difference 

A = £l -e (12) 

between the lowest two levels is not sensitive to the vari- 
ation of (3l- In Fig. 4, we show the energy separation of 
the two lowest levels, A, as a function of (3^. We find the 
interesting result that A(/?^) is almost flat at / = 0.5 
(0.011 < A(p L )/Ej < 0.0135) when < (3 L < 0.85. 
These features indicate that, even with a large loop in- 
ductance of (3l = 1, in the vicinity of / = 0.5 the two 
lowest eigenstates (denoted by |0) and |1) for the ground 
and the first excited states, respectively) remain suitable 
basis states for a flux qubit. Within the subspace of qubit 
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FIG. 4: Energy difference A between qubit states |1) and 
|0) as a function of /3l for / = 0.5. Here A is in units of 
Ej. Notice that the energy difference varies ~ 0.003-Ej when 
varying the loop inductance L. 

states spanned by |0) and |1), the Hamiltonian is reduced 
to 

fr = ei|l)<l|+e |0)<0|. (13) 

If the average energy (e\ + £o)/2 is chosen to be the new 
zero-point energy of the flux qubit, the Hamiltonian can 
be further expressed as 

H=\a Pz , (14) 

where p, = |1)(1| - |0)(0|. 

C. Comparision with other works 



the external flux <E> e is applied within the loop A\L C B\A\. 
Also, the circuit is designed in such a way that the mutual 
inductance between loops A\L C B\A\ and A\B\B 2 A 2 A\ 
may be ignored. This is achieved when only a small frac- 
tion of the flux generated by one loop passes through 
the other. (If this were not to be the case, the interbit 
coupling can still be achieved by the common inductance 
L c , but the interaction Hamiltonian takes a more com- 
plicated form.) Phase drops through the three Josephson 
junctions of the ith flux qubit are constrained by 

<f>f - <t>f + 4>f + 2tt [/ + {UU + IjL c )/<f> ] = 0, (15) 
where i,j = 1,2 (i ^ j), and 

L\ = L c + L/,1, 

L 2 = L c + L 12 + L b2 . (16) 
The total supercurrent through L c is 

I = h+h, (17) 

where 

Ii = lot sin(0 pi + (f> mi ), (18) 
with I 0i = and 

<t> mi = \{<t>¥ -<$)■ (19) 

The Hamiltonian of the two coupled flux qubits can be 
written as 



In Ref. |24j, the effects of the loop inductance in a 
flux qubit are considered using a perturbation approach, 
where the Hamiltonian is expanded into three parts: an 
inductance-free Hamiltonian, an inductance-related har- 
monic oscillator term, and a small correction term. This 
perturbation method is valid for /3l -C 1 since the cor- 
rection term is proportional to the loop inductance of the 
flux qubit. Instead of using the perturbation approach, 
we numerically solve Eq. (11) to obtain the eigenvalues 
and eigenstates of the system. This numerical method 
allows us to extend our study to the regime of 0l ~ 1, 
where the lowest two eigenstates of the system can still 
be used for achieving a qubit. Using the experimental 
valued I ~ 0.5 p,A, this regime corresponds to a loop 
inductance of L ~ 1 nH. 



H = H 1 +H 2 + H I . (20) 

Here Hi is the Hamiltonian of the ith isolated flux qubit, 
with loop inductance Li and circulating supercurrent Ii, 
which has the form in Eq. (5) but with /' replaced by 

./•: ./• • (2i) 

Also, Ej, C and a are replaced by Ej, and a^. The 
interaction Hamiltonian is 

2 

H^Lchh-Y^aiEfu,, (22) 
i=i 

where 



III. COUPLED FLUX QUBITS n » = cos (7> + 2TrIjL c /%) - cos 7i , (23) 

with 

7l = 2tt/; + 20 im , (24) 



A. The model 



To couple two flux qubits, we use a common induc- 
tance L c shared by these two qubits [see Fig. 1(b)]. Here and i ^ j. 
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When 0Li = 2irIoiL c /$o <C 1, Hj is approximated by 
Hi = -L c hh (25) 

because 



(26) 



in this case (see the Appendix). Within the qubit-state 
subspace of the ith isolated flux qubit, Hi is reduced to 



where 



H = ±A lP <i\ 



p« = |l i )(l i |-|0,)(O i |. 



(27) 



(28) 



In the vicinity of / = 0.5, because the supercurrents Jj at 
states and |0j) have equal magnitudes but opposite 
directions, Ii can be written as 



Ii=a iP ^ +b i \l i )(0 i \ + b*\0 i )(l i 



(29) 



LU 



LU 



where 




0.51 



0.51 



k = (U\h\<k). (30) 

Because the supercurrent Ii at state (i.e., a,) is pro- 
portional to the slope of the energy level that corresponds 
to state \li) with respect to / (see, e.g., Ref . il9f) . it falls 
to zero at the symmetric point, / = 0.5, where the level 
becomes flat. Also, our numerical results show that bi 
becomes a real number at / = 0.5. Thus, we can rewrite 
Ii at f = 0.5 as 

Ii = bi P ®, (31) 

with 

pW=|l i )<0 i | + |Q < )(l i |. (32) 

For [3 n <C 1, i.e., the common inductance L c is very 
small, the Hamiltonian at / = 0.5 can be cast to 

n = j^\*i&-x&V?>, (33) 
«=i 

with 

X = L c btb 2 . (34) 
It is clear that the interbit coupling persists at / = 0.5. 

B. Energy spectrum 

Figure 5 shows the energy spectrum of the two coupled 
flux qubits around / = 0.5. In order to realize fast two- 
bit operations while keeping the leakage from the qubit 



FIG. 5: (Color online) Energy levels of two coupled flux qubits 
versus reduced flux / for Lt\/L c = 0.1 and (L12 + L^)/ L c = 
0.2. The parameters Pli = 2nIoiL c /$o are (a) (3li = Pl2 = 
0.03, (b) f3 L1 = 0.03, I3 L2 = 0.04, (c) /3 L1 = f3 L2 = 0.07, and 
(d) Pli = 0.07, Pl2 = 0.1. Here the energy E is in units of 
-Ej 1 '. Near / = 0.5, the energy levels are robust with respect 
to the asymmetry between /3li and /Mis- 



states to other higher energy states small, we choose the 
interbit coupling strength to be comparable to the energy 
difference, at / = 0.5, between the basis states and 
|0j) of each qubit. As shown in Figs. 5(a) and 5(b), 
the energy spectrum remains similar in the vicinity of 
/ = 0.5 when the two flux qubits have different values of 
parameters. Furthermore, the two higher energy levels, 
€3 and £4, in the first four energy levels (i.e., with 
k = 1 to 4) of the two coupled flux qubits are flat in a 
relatively broad range around / = 0.5; this flat region 
is much broader than the corresponding flat-energy-level 
range of the single flux qubit around / = 0.5. The flux- 
independent level £3 in Fig. 5(a) corresponds to a singlet 
eigenstate, while other three levels correspond to triplet 
eigenstates. As expected, the transitions between this 
singlet state and other three triplet states are not allowed 
by the microwave perturbation [cf. Fig. 6(a)]. 

When the interbit coupling increases further, the flat 
region for both levels £3 and £4 widens for two qubits 
having identical parameters [see Fig. 5(c)], but £3 and £4 
become much different in this region when the two qubits 
are not identical [see Fig. 5(d)]. Moreover, it can be seen 
that, at / = 0.5, the gap between levels £1 and £2 and 
that between £3 and £4 become narrow when increasing 
the interbit coupling. 
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At / = 0.5, the first four energy levels, e&, k = 1 to 4, 
of the coupled flux qubits can be approximated by 

?2 = ~\ E B, e 4 = -E A , (35) 

where 

S A = [(A 1 + A 2 ) 2 + 4x 2 ] 1/2 , 

i? B = [(A 1 -A 2 ) 2 +4x 2 ] 1/2 . (36) 

The gap between levels e 2 and £3 is -Eb, which increases 
with x- The gap between levels t\ and e 2 and that be- 
tween £3 and 64 are given by Ea — -Eg. Figures 5(a) 
and 5(b) correspond to x ~ Ai; in Fig. 5(a) where 
Ai = A 2 = A, the two equal gaps, Ea — Eb, at / = 0.5 
are (\/2 - 1)A. When X further increases, the value of 
Ea — Eb decreases; namely, the two equal gaps become 
narrow [cf. Figs. 5(c) and 5(d)]. 

In the case of Fig. 5(a), because 2-Kl 0i L c /§ = 0.03, 
the common inductance is 

L c ps 20 pH 

if the critical currents igi are equal to the experimen- 
tal valued Iq ~ 0.5 /iA. Such a small inductance is ex- 
perimentally realizable, e.g., using a loop of diameter 
d ~ 16 fjxa. Also, our numerical calculations show that 
hi ps O.66/0 at / = 0.5. The interbit coupling is thus of 
the order 

X = LMhi ~ OMSEj, (37) 

which is equal to the energy difference A at / = 0.5 
of the single flux qubit with /3l — 0.03. The typical 
two-bit operation time related with the interbit-coupling 
strength is 

t 2 ~ h/ x . (38) 

Because x ~ A when L c ps 20 pH, 

t 2 ~ n = K/A, (39) 

where t± is the one-bit operation time. Therefore, the 
corresponding two-bit operation is as fast as the one-bit 
operation. 

C. Comparison with other works 

Spectral results similar to the ones shown in Figs. 5(a) 
and 5(b) were also obtained by Storcz and Wilhelm 21 
and by Majer et alml using simpler model Hamiltonians 
for two coupled flux qubits. However, because a different 
setup is used in Ref. |22 for coupling the two qubits, the 
four energy levels are flipped as compared to ours. 

In our proposed setup, the shared part of the loop is 
on the same side of each qubit and, as given in Eq. (37), 



X > 0. Since the negative coupling term in Eq. (33) fa- 
vors the parallel arrangement of the (pseudo)spins, the 
interbit coupling here is ferromagnetic. In contrast, in 
the setup of Ref. |22j, the shared loop is on the opposite 
side of each qubit. This gives rise to a negative coupling 
parameter x- Therefore, the interbit coupling there is 
antiferromagnetic because the (pseudo)spins tend to ar- 
range antiparallel to each other. 

Here we use the common loop inductance to couple 
two flux qubits, while different setups are proposed in 
Refs. and I20I where the mutual inductance is used 
for coupling flux qubits. Since the mutual inductance 
is always smaller (sometimes can be much smaller) than 
the loop inductance of each flux qubit, if the two flux 
qubits are coupled via the mutual inductance (instead of 
the common loop inductance in our approach), a larger 
(sometimes much larger) loop is required for each qubit 
to produce a strong enough interbit coupling. There- 
fore, due to the larger loop in each flux qubit, the system 
would experience more serious flux noise. This is a very 
significant difference between these approaches. 

In our approach, we treat both single and coupled flux 
qubits using more realistic models including the loop in- 
ductance. If our more general theory were used to deal 
with the different setup proposed in, e.g;, Ref. 12^, the 
simple model Hamiltonian used in Ref. |22j would be an 
approximation of the more general theory. The simpler 
model can explain the current experiment in Ref. l22l 
and so does the general theory. Moreover, the more gen- 
eral theory presented here can explain additional features 
that would be relevant for future experiments, while the 
simple model might not. For instance, even for the sin- 
gle flux qubit, the simple model cannot tell how large the 
loop should be to break down the qubit. Moreover, the 
simple model involves only two levels for each qubit. It 
cannot explain phenomena related with the state transi- 
tions from these two levels to higher ones. However, the 
general theory can do it. 



IV. STATE TRANSITIONS 

A. Transition matrix elements for single and 
coupled two qubits 

When a microwave field with an appropriate frequency 
uj is applied through the superconducting loop of the sin- 
gle flux qubit, a transition between two states occurs. 
Now, the total flux within the loop is <fr e + where 

$ f (t) = $ x cos(wt + 6>) (40) 

is the microwavc-ficld-induced flux through the loop. For 
a weak microwave field, the single flux qubit experiences 
a time-dependent perturbation 

H'(t) = -I$x cos(ujt + 9), (41) 



7 




0.49 0.50 0.51 0.49 0.50 0.51 



f f 



|0) — > |2) is forbidden. Also, it can be seen that | ioi I is 
not sensitive to the variation of /3l, while | io2 1 and \ti2\ 
are slightly reduced when increasing (3l- This observa- 
tion is consistent with the energy spectrum in Fig. 3, 
where the gap between the lowest two levels and 1 is 
not significantly changed, but the gap between levels 1 
and 2 slightly increases with (3l- In Figs. 6(b) and 6(c), 
we show the flux dependence of \tij\ for all possible tran- 
sitions in the coupled flux qubits. When the two flux 
qubits have the same parameters, the transitions 

lei) - |ea), 

|e 2 ) - lea), (44) 
lea) - |c 4 ) 

are forbidden because \tij\ = [see Fig. 6(b)]. However, 
they are allowed (except for some specific values of /) 
when the parameters of the two flux qubits are different 
[see Fig. 6(c)]. These properties are attributed to the 
changes in the symmetry of the qubit states . In particu- 
lar, when / = 0.5 , \ti2\ has the largest value, while |t24| 
has a smaller value and others are either zero or much 
smaller. 



B. One- and two-bit operations implemented via 
microwave fields 



FIG. 6: (Color online) (a) Moduli of the transition matrix 
elements tij between single-qubit states \i) and \j) versus re- 
duced flux /. (b) and (c) Moduli of the transition matrix 
elements between coupled-qubit states |e») and \ej) versus / 
for L b i/L c = 0.1 and (-L12 + L b2 )/L c = 0.2. Here \Uj\ is in 
units of Iq<&x in (a) and Joi$x in (b) and (c). Notice that, in 
(b) and (c), some transitions for the coupled two qubits are 
sensitive with respect to the asymmetry of /3li and /3l2- 

and the transition matrix element iy between states \i) 
and \j) is given by 

Ui = (i\I<S>x\j). (42) 

Similarly, when the microwave field is applied through 
the left loop A\L C B\A\ of the coupled flux qubits, the 
transition matrix element Uj between the coupled-qubit 
states \ei) and \ej) is then 

= <ei|7$x|ei>. (43) 

Notice the difference between Eqs. (42) and (43), though 
they look similar to each other in expression; in Eq. (42), 
\i) is an eigenstate of the single flux qubit and I is the 
circulating current in the qubit loop, while / in Eq. (43) 
is the total current I = I\ +I2 in the coupled-qubit circuit 
and |ej) is an eigenstate of the coupled two flux qubits. 

Figure 6(a) presents the flux dependence of for 
transitions |0) -> |I>, |0) -> |2), and |1) -> |2) in a sin- 
gle flux qubit. Because of the symmetry of the wave 
functions, |io2 j = at / = 0.5, and thus the transition 



For a single flux qubit with (3l = 0.03, the energy 
difference A between states |1) and |0) is 0.0I29I-E/ at 
/ = 0.5. Using an experimental valued for the critical 
current Iq ~ 0.5 //A, we obtain Ej ~ 1.03 meV. The 
energy difference of 0.0I291i?j corresponds to a gap of 
v w 3.2 GHz. The one-bit operation can be implemented 
using a resonant microwave field. For a weak driving 
field, the Rabi frequency fioi is given by \toi\/H- The 
typical switching time is isw — Tr/^oi when the states 
|0) and |1) flip. For instance, because |ioi| ~ 0.66Jq$x 
at / = 0.5, the switching time isw is about 3 ns for 
I $x ~ I jueV. If the leakage from these two states to 
others is small, one can realize a fast one-bit operation, 
e.g., with a switching time 

^sw = 7^— ~ 10v (w 3 ns), 

"01 

by increasing the microwave-field intensity. 

Let the energy difference between states |0) (|1)) and 
|1) (|2)) be huiQi (I1UJ12). When the field is tuned to be 
resonant with the transition |0) — > |1), the ratio of the 
transition probabilities between |1) — > |2) and |0) — > |1) 
can be estimated as 

P12 _(^i2\ 2 sin 2 (Or/2) 
P01 V / sin 2 (i2 01 r/2)' 

where 

n = {n 2 12 + (u Jl2 ~Lu 01 ) 2 } 1 / 2 , (46) 
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^12 = |ii2|/^, and r is the duration of the microwave- 
field pulse. When 



7T/U 



01 



I0u~ 



using the numerical results Huqi = 0.01291-E,/, hu>i2 
0.18763.BJ, and |*i 2 /*ox] » 0-38 at / = 0.5, we have 



P12 
Pox 



1.5 x 10" 



(47) 



This implies that the leakage to other states is small for a 
fast one-bit operation implemented via a microwave field. 

Corresponding to Fig. 5(a), |ei) and |e 2 ) at / = 0.5 are 
approximated by 



1 



^(|oi) + |io», 



(r, |Q0) + |11)), 



with 



V = 



A + (A 2 + x 2 ) 
X 



2^ 1/2 



(48) 



(49) 



Initially preparing the system at the (entangled) ground 
state |ei), one can produce the maximally entangled state 
|e 2 ) using a microwave-field pulse of duration r — 7r/r2 12 , 
where the Rabi frequency f2i2 is given by |ii2|/?i for a 
weak driving field. 

At / = 0.5, we have hu>i 2 = 0. 00528£ , j, hu 2 4, = 
0.03124/5./, and |*24/*12| « 0.41. When the microwave 
field is in resonance with the transition |ei) — > leg) at 
/ = 0.5, 



for 



« 4.4 x 10- 6 

P12 



T = 7r/f2l2 ~ 207r/wi 2 . 



(50) 



Because the state leagage is very small, a fast two-bit 
operation can also be implemented using a microwave 
field. 



SUPERCURRENTS AND QUANTUM 
MEASUREMENT 



The circulating supercurrents flowing through the in- 
ductance L or L c are different for different eigenstates. 
This property can be used for implementing a readout of 
the qubit states. For a single flux qubit, around / = 0.5, 
the supercurrents I at eigenstates |0) and |1) (i.e., (0|/|0) 
and (1|/|1)) have equal magnitudes but opposite direc- 
tions [see Fig. 7(a)]. During quantum measurement, one 
can switch on the flux transformer to couple the induc- 
tance L with a dc-SQUID magnetometer [cf. Fig. 1(a)] 
to distinguish the two eigenstates of the qubit because at 
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FIG. 7: (Color online) (a) Supercurrents I versus reduced flux 
/ at eigenstates |0) and |1) of a single flux qubit for /3l = 0.03. 
(b) Supercurrents I\ versus / at eigenstates \eh), k = 1 to 
4, of two coupled flux qubits for the symmetric circuit with 
Pli = Pl2 = 0.03. (c) Supercurrents 7i, I2, and I versus / at 
eigenstates \ek}, k = 1 to 4, of two coupled flux qubits for an 
asymmetric circuit with /3li = 0.03 and 0L2 = 0.04. Here we 
choose Lb\/L c = 0.1 and (L12 +iw)/ic = 0.2 for the coupled 
flux qubits. The supercurrents are in units of 7o in (a) and J01 
in (b) and (c). For the coupled two qubits, the total current 
I = Ii + I2 is robust with respect to the asymmetry of (3li 
and Pl2- 



these two states the supercurrents / through L generate 
two different fluxes in the SQUID loop of the magnetome- 
ter. In general, if the single flux qubit is at the superpo- 
sition state ci 1 1) +co|0), the measurement will show that 
the qubit has probability \a\ 2 at the eigenstate \i), where 
i = 0,1. For the two coupled flux qubits, the supercur- 
rents through the common inductance L c take different 
values at its four eigenstates. 

Similar to the single flux qubit, a switchable flux trans- 
former can be used to couple L c and the SQUID loop 
of the magnetometer for reading out the coupled-qubit 
states because the supercurrents / at different eigenstates 
contribute different fluxes in the SQUID loop of the mag- 
netometer. The supercurrents I\ at the four eigenstates 
of the coupled qubits are shown in Fig. 7(b) for two flux 
qubits having identical parameters. Since I\ = I 2 in this 
case, the total supercurrent I is 2I\. When the param- 
eters of the two flux qubits become different, the total 
supercurrents I look similar to those in Fig. 7(b), but I\ 
and I 2 (which flow through the Josephson junctions of 
the qubits) change drastically [cf. Fig. 7(c)]. Also, it can 
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be seen that at the eigenstates of the system the circulat- 
ing supercurrents in both single and coupled flux qubits 
fall to zero at / = 0.5. To read out the qubit states, one 
can shift the system away from this point. 



of the two qubits change from being initially equal to 
each other and then making these different. 

Acknowledgments 



VI. DISCUSSION AND CONCLUSION 

For the charge qubits coupled by LC-oscillator 
modest or by an inductance^ the inductances pro- 
posed to be used are ~ 3.6 /zH or ~ 30 nH, respec- 
tively, for a two-bit operation ten times slower than the 
typical one-bit operation. An inductance for coupling 
charge qubits similar to that in Ref . Il3l particularly, has 
a value much larger than L c (« 20 pH) for coupling flux 
qubits. ft is difficult to fabricate in a small size with- 
out introducing a strong coupling with the environment. 
Because two-bit operations are much slower than one-bit 
ones in the inductively coupled charge qubits, an effi- 
cient scheme is thus required to minimize the number of 
two-bit (as opposed to one-bit) operations to obtain a 
conditional gate, 1:4 However, for inductively coupled flux 
qubits, the above limitation in using two-bit operations 
for constructing a conditional gate is removed since two- 
bit operations can be as fast as one-bit ones. In this 
case, any schemes for constructing conditional gates be- 
come efficient by minimizing the number of operations 
that are used (either one- or two-bit). Note that the 
common inductance of L c w 20 pH can produce a strong 
interbit coupling. As a result, two-bit operations as fast 
as one-bit ones can be achieved. This common induc- 
tance is comparable to the loop inductance, L ~ 10 pH, 
of the single flux qubit currently realized in experiments. 

To couple several flux qubits, the inductances of all 
loops involved could be small, comparable to the loop 
inductance of a single flux qubit currently realized in ex- 
periments. This is the case we studied in the present 
paper, where two coupled flux qubits are considered. If 
a number of flux qubits are coupled, the inductances of 
some loops will become larger, but the common or shared 
inductance for producing the interbit coupling can still be 
chosen small (about 20 pH). If the circuits except for the 
line A\L c Bi (corresponding to the common or shared in- 
ductance) could be screened from the environment (that 
is a big challenge for experimentalists for sure) , the main 
noise would be due to the small common or shared in- 
ductance. 

In conclusion, we have proposed an experimentally re- 
alizable method for inductively coupled flux qubits that 
can achieve two-bit operations performing as fast as one- 
bit ones. We treat both single and coupled flux qubits 
with more realistic models including the loop inductance. 
Moreover, we show that the coupled flux qubits have 
novel flux-dependent behaviors in the transitions between 
states. We find that the forbidden transitions in the cou- 
pled two flux qubits become allowed (except for some 
specific values of the external flux) when the parameters 
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APPENDIX: SERIES EXPANSION OF THE 
INTERACTION HAMILTONIAN 

The supercurrent Ii flows through each of the three 
Josephson junctions in the zth flux qubit, so Jj can also 
be written as 

Ii = ail 0i sin (j>$ 

= -aj 0l sin {2tt(/; + IjL c /<S> ) + 20 mi } , (A.l) 

where i,j = 1,2 and i ^ j. Taking advantage of this re- 
lation for Ii , one can expand the interaction Hamiltonian 
(22) as 



H I = -\L c I 1 I 2 -Y J ^Efi 



(A.2) 



i=l 



where 
A= 1 



E 

i=l 



-A 



Li 



1^ 

hi 



15 



0h 



loi 



and 



(A.3) 



& = - cos(27t/^ + 2</> lm ) 



Li 



with 



P Ll = 2nI 0i L c /<$>u < tt/2. 



(A.4) 



(A.5) 



The term —\L c I\l2 in Hi produces an interbit coupling 
between flux qubits 1 and 2, while onEj^i slightly mod- 
ifies the energy levels of the ith flux qubit. 
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